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ABSTRACT 

We have found an optical/X-ray counterpart candidate for the bright, but presently unidentified, Fermi source 
2FGL J131 1.7— 3429. This counterpart undergoes large amplitude quasi-sinusoidal optical modulation with a 
1 .56h (5626s) period. The modulated flux is blue at peak, with T e g « 14,000 K, and redder at minimum. Super- 
imposed on this variation are dramatic optical flares. Archival X-ray data suggest modest binary modulation, 
but no eclipse. With the 7-ray properties, this appears to be another black-widow-type millisecond pulsar. If 
confirmation pulses can be found in the GeV data, this binary will have the shortest orbital period of any known 
spin-powered pulsar. The flares may be magnetic events on the rapidly rotating companion or shocks in the 
companion-stripping wind. While this may be a radio-quiet millisecond pulsar, we show that such objects are 
a small subset of the 7-ray pulsar population. 
Subject headings: gamma rays: stars — pulsars: general 



1. INTRODUCTION 

In the most recent Fermi Large Area Telescope (LAT) cat- 
alog of 1873 0.1 — lOOGeV sources, over 1170 have sta- 
tistically reliable c ounterpart identifications at lower energy 
dNolan et alj 120121) . The bulk of these identifications are 
blazars and spin-powered pulsars, with young pulsars (both 
radio-selected and Geminga-like 7-ray selected) dominating 
at low Galactic latitude. An additional population of mil- 
lisecond pulsars extends to high latitude among the blazar- 
associated sources. Identification progress is especially im- 
pressive for the bright, well localized sources that have been 
studied since early in the Fermi mission. Here we update 
lAbdo et ail d2009al) . selecting 'bright' sources from the 2FGL 
catalog as those with > 20cr detection significance and time- 
averaged energy flux Fe > 3 x 10~ n ergcm 2 s _1 - there 
are 249 such sources. All but six presently have lower en- 
ergy identifications: blazars, pulsars and a few binaries. The 
handful of sources not yet associated with one of these source 
classes provide the best prospect for new types of 7-ray 
emitter. We have initiated a campaign to characterize these 
unidentified sources. The first result from this effort is the dis- 
covery of dramatic optical and X-ray variability for the high 
latitude = 62°) source J2339-0533, implying that it is 
a short period 'black- widow' - t ype binary millise cond pulsar 
(MSP) ((Romani & Shawll20miKorig. et al.ll2012h . 

Here we report progress on the next high latitude unidenti- 
fied source in this set, 2FGL J 1 3 1 1 .7-3429 (hereafter Jl 3 1 1 ) 
at |6| = 28°. With a 43er detection significance (the high- 
est among the 2FGL unidentified sources) and an energy flux 
of -Fo.i-iooGeV = 6.2 x 10 -11 ergcm~ 2 s _1 (the second 
brightest UnlDed 2FGL), this is a top candidate for follow- 
up. The 'Variability index' value 19.5 and 'Curvature Sig- 
nificance'=6.3 mark this as a steady source with a substan- 
tial spectral cut-off: a prime pulsar can didate. It has been 
searched for < 30 Hz 7-ray pulsations (eg. lAbdo et alJ2 009b) 
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and for radio pulses down to ~ ms periods dRav et alj|2012l) . 
with no detection. Thus it is unlikely to be an isolated young 
pulsar or a radio-loud MSR We describe here an optical cam- 
paign to identify and characterize a counterpart. 

2. OBSERVATIONS 

Initial exposures of the J131 1 error region were taken with 
the MiniMo camera at the 3.6 m WIYN telescope on Feb. 17- 
18, 2012. Exposures were 180 s using the Gunn filter set (2 x 
g, r, i on 2/17; 5 x g, r on 2/18). Observations were perforce 
at airmass 2.5 — 3, seeing was variable and conditions were 
non-photometric. 

A more extensive image sequence was taken with the 
SOAR Optical Imager (SOI) at the 4.2 m SOAR telescope on 
March 21-23, 2012. Here 120 s exposures using SDSS filters 
were taken, covering one color each night (61 x g' on 3/21, 
74 x r' on 3/22 and 57 x i' on 3/23). The camera was binned 
2x2, allowing fast lis read-out and high observation effi- 
ciency. Again image quality was highly variable (partly due 
to loss of wavefront sensor control of the primary). 

We searched image archives for exposures covering 
J1311, finding a sequence (program 086.D-0388) of 170s 
VLT/VIMOS exposures taken on March 4 (2 x [3B, 3V, 31}) 
and March 11, 2011 (3 x [3B, 3V, 31}). VIMOS is an array 
camera with 2' gaps between the four CCD chips, so the five 
pointings were offset to cover the full J131 1 error region. 

In examining the images for variable objects, we paid 
particular attention to sources coincident with detections in 
archival X-ray data (see below). One such star was blue and 
highly variable on short timescales. The NOMAD position of 
this likely counterpart is 13 1 1 45.741, -34 30 29.96 (2000.0). 
While we ensured that this star was covered in all MiniMo 
and SOI exposures, we found that the counterpart was in the 
array gap for the first VIMOS pointing on March 3, 201 1 and 
at the extreme edge of a chip, suffering major vignetting, in 
the second pointing on March 1 1, 201 1. 

To quantify the variations, we performed simple aperture 
photometry of the counterpart and surrounding field stars in 
all frames. Since the bulk of the imaging was in the SDSS 
filter set, we attempted to reference all photometry to this 
scale. The MiniMo data were calibrated using gri images 
in SDSS fields at low airmass; we transferred this calibra- 



2 



Romani 



50 




Binary Phase 



FIG. 1.— SOARAVIYN/VLT photometry of J131 1.7— 3429. The abscissa 
is the phase relative to phase (superior conjunction, secondary flux maxi- 
mum) just before the first observation of each plotted night, using the best-fit 
Pl> =5626.00 s period. Symbols and colors indicate the various telescopes 
and observing bands. Photometry from the SOAR i' flare is connected with 
a solid line; points for each VLT night are connected with a dash line. 

tion to stars in the J 1 3 1 1 field. Unfortunately, the non- 
photometric WIYN conditions and pass-band differences be- 
tween the SDSS and Gunn sets make the substantial atmo- 
spheric extinction corrections for J 1 3 1 1 uncertain. This is 
not a severe problem for the field stars, which were mostly 
late type with g' — r' = 0.7 — 1.5. However the target's 
color is as blue as g' — r' = —0.5 and varies appreciably 
with magnitude. Thus while the relative photometry for each 
color sequence has very small statistical errors, there may be 
a systematic offset between the WIYN and SOAR points as 
large as 0.2 mag, especially in g' . For the VIMOS images, 
we used NOMAD field stars to establish zero points in B and 
V. NOMAD-detected stars were generally saturated in the VI- 
MOS I images, so we established the I scale using stars from 
the DENIS catalog. 

We show the photometry for each night in Figure 1, con- 
verting to fluxes and plotting against estimated binary phase. 
The SOI g', r' data show dramatic quasi-sinusoidal variabil- 
ity by a factor ~ 30 x with a period of ~ 1.5 h. In i' an 
initial flux minimum was followed by a large ~ 5x flare, 
with rapid fluctuations during return to quiescence. The g', r' 
fluxes also show variability, especially near maximum. The 
WIYN data are sparse and have larger statistical errors, but 
do follow the periodic modulation. Most WIYN fluxes ap- 
pear 10-20% lower than SOAR at the corresponding phase, 
suggesting an underestimate of ~ 0.1 mag for the extinction 
correction. However, given the apparent non-periodic vari- 
ability this cannot be firmly concluded. The VLT fluxes from 
March 1 1 pass through two minima; those from March 4 show 
one minimum, but appear to have ~ 8fi]y additional flux. 

2.1. Orbital Period Estimate 

The intermittent flaring activity of J131 1 compromises sim- 
ple periodogram or Fourier methods for refining the period. 



However, the combination of SOAR, WIYN and VLT expo- 
sures give photometry information on 1-3, 7, ^35, and ^350 
day time scales, which suffices to resolve possible aliases. Af- 
ter referencing all observation mid-point times to the solar 
system barycenter, we can constrain the phase of the 2012 
epoch minima to 8(f) « 0.02, giving a best fit period of 
Pb = 5626.00 ± 0.02s. The epoch of (quiescent) maximum 
light associated with our SOAR r' measurements (presum- 
ably pulsar superior conjunction, see below) is barycentric 
MJD 56009.1795 ±0.0013. 

The approximately sinusoidal nature of the light curve and 
very large modulation amplitude can be best interpreted in 
terms of strong companion heating. The inferred heating lu- 
minosity is appreciably above any X-ray flux (see below); we 
infer that accretion power is not dominant and that the heat- 
ing source is invisible in the optical and faint in the X-ray. 
We conclude that this is a 'black-widow' type pulsar binary 
and is the likely counterpart of the 7-ray source. In this it 
is remarkably similar to the black widow candidate counter- 
part for J2339— 0533. But the ~ 3x shorter orbital period is 
unprecedented; this is an extreme black widow system. 

2.2. Archival X-ray Light Curve 

2FGL J13 11.7-3429 has been observed by both the CXO 
and Suzaku satellites. One exposure from each facility is 
available in the NASA archives, additional exposure has been 
taken but is not publicly available. The counterpart was de- 
tected in both data sets and well localized by CXO. We ana- 
lyzed both observations for X-ray spectrum and variability. 

The Suzaku observation (ObsID 804018010; Kataoka, PI) 
started on August 04 2009 (MJD 55047.2290) and had 33.4ks 
of on-source time over 95 ks (~ 17 binary orbits) with ex- 
posure in XIS0,1,3. A comparably bright source lies ~ 1.6' 
from the target, so we extracted source counts from a 1' radius 
aperture. Figure 2 (left) shows the combined XIS count rate, 
both direct aperture counts (histogram) and the background- 
subtracted aperture-corrected count rate (points). The vertical 
dotted lines mark the phases of maximum optical light; we 
see that orbital modulation is difficult to measure since Pi, is 
so close to that of satellites in low Earth orbit. A clear X-ray 
flare is seen 5000 s after the observation start. 

An ACIS-I exposure of 19.8 ks live-time (ObsID 11790; 
Cheung, PI) was taken starting MJD 5527.6652. Unfortu- 
nately, as for VLT observation 1, the counterpart was directly 
in the gap between ACIS-I chips. However, the CXO dither 
did provide some (^25% effective live-time) exposure and a 
source coincident with the optical counterpart was well de- 
tected with 60 counts. There is no strong evidence for secu- 
lar X-ray variability on 3000 s timescales. The exposure only 
covers ~ 3.5 orbits. Figure 2 (right) shows the folded X-ray 
light curves from the Suzaku and CXO data sets. After the 
early flare in the Suzaku data is excluded there is only weak 
evidence for orbital modulation; note there is an exposure gap 
for phases 4>b = 0.7 — 0.8. Similarly the CXO data show no 
strong peak or eclipse. In both light curves there seems to be 
a slight excess at phases 0.7-1.1. This may indicate additional 
weaker flaring activity near superior conjunction. 

Only crude spectral constraints can be obtained with the 
limited counts, especially after selecting for Suzaku photons 
outside of the flare period. If we fit an absorbed power 
law to the X-ray data we find that Nh is essentially uncon- 
strained (0.5 ± 6 x 10 21 cm -2 ) but consistent with the lim- 
ited extinction in this high \b\ direction. We therefore fix 
at the Nh = 5 x 10 20 cm~ 2 inferred from the dust maps 
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FIG. 2. — Left: Raw aperture count rate (histogram) and background- and aperture- corrected counts (points) during the Suzaku exposure. Times of optical flux 
maximum are marked. Right: Two periods of the X-ray orbital light curve of J1311.7— 3429. Statistical errors are shown for the exposure-corrected bin count 
rate. For Suzaku, we show the folded data both with and without the bright flare. Horizontal error bars show the extrapolated phase uncertainty. 
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FIG. 3.— Folded SOAR g'r'i' light curves of J131 1.7-3429 (two periods). 
Statistical errors are shown during the first period only and the large i' flare is 
plotted during the second period only. The curves represent a simple heated 
secondary model from ELC; a second copy of the i' curve with a constant 
flux offset is also shown for comparison. 

(HEASARC NH tool), and fit for the power law index. Us- 
ing CIAO the CXO data fit to T = 1.3 ± 0.5. XSPEC 
fits to the Suzaku XIS0,3 give T = 1.6 ± 0.8. Fixing at 
T = 1.5, we obtain a quiescent unabsorbed flux for the source 
of /o.3-iofceV = 2.4 ± 0.6 x 10~ 13 ergcm -2 s _1 . A deeper 
exposure, allowing a better light curve and a orbit-dependent 
spectral study is clearly of interest. 

3. SYSTEM MODELING 

Figure 3 shows two periods of the SOI light curves folded 
on the ephemeris of §2.1. The large quasi-sinusoidal modu- 
lation shows that the light curve is dominated by heating of 



the secondary. The typical g — r w —0.4 near = sug- 
gests T cff » 14, 000 K. By phase = 0.35 the color has 
reddened to g - r > -0.1 (T eS » 8000 K). At flux mini- 
mum (<fi ~ 0.5) large (up to Am = 3) fluctuations dominate. 
Since relatively small fluctuations can overwhelm the quies- 
cent flux, these prevent any meaningful measurement of color 
at flux minimum. 

We have attempted to model the 'quiescent' light curves 
with the ELC code (lOrosz & Hauschildtl [20001) . A formal 
fit is not indicated as there are many positive 'flare' excur- 
sions from any reasonable light curve. However compari- 
son with a range of models shows that a few of the param- 
eters are significantly constrained. For example, matching 
the characteristic g' — r' at maximum requires heating by 
a primary luminosity (modeled here as an isotropic 'X-ray' 
flux) of Log[L x (erg/s)] = 36.0 ± 0.3. Also, to obtain ad- 
equate modulation while avoiding an X-ray eclipse, requires 
60° < i < 80°. At this point other parameters are poorly 
constrained, although large mass ratios and a near-Roche- 
lobe filling secondary are preferred. However, the basic ELC 
model seems inadequate for a good fit; a particular challenge 
is to reproduce the relatively fiat maxima and broad 'shoul- 
ders' in the g' and r' curves at <fi ~ 0.25, 0.75. This may be 
due to limitations of the atmosphere tables. One can substan- 
tially improve the fit by modifying the limb-darkening law, 
but for blackbody surface spectra broad light curves (with ad- 
equate heating) require unphysical negative limb darkening 
coefficients. Another, more attractive, possibility is that the 
distribution of surface heating differs from the isotropic form 
assumed here; the pulsar may primarily heat a portion of the 
companion star. An illustrative set of curves is shown in Fig- 
ure 3, assuming a mass ratio q = 50, orbital semi-major axis 
a = 0.9-R©, inclination i = 65°, heating flux 10 36 ergs _1 and 
a T e ff = 4000 K companion reaching 0.99 of its Roche lobe 
radius. The observed i' flux evidently has a substantial added 
flux even before the bright flare, at the epoch observed. 

4. DISCUSSION 

At this short orbital period, the mean density in the sec- 
ondary Roche lobe is (p2) = 46gcm -3 , so the secondary 
could be a main sequence star of < 0.15M Q . To reach 
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this period, however, the system very likely passed through a 
common envelope phase, and would have a helium-rich core. 
The strong heating and near-Roche lobe filling suggest a sub- 
stellar secondary, but additional observations are needed to 
determine the state of the secondary core and envelope. 

A very rough estimate of the source distance can be made 
from the g' — r' = —0.4 color and g' = 20.2 flux at max- 
imum, if we assume that the star nearly fills its Roche lobe 
at pa 0A6a(M 2 /M Tot ) 1 / 3 pa 0.1R e . Such color indicates 
T c s ~ 14,000 K (type ^B6), corresponding to main sequence 
M g = —0.43 and R pa 3.0 Rq. Correcting for the effective 
area and the small (A g = 0.25) Galactic extinction in this 
direction, we get d pa 3.9(<7/50) _1 / 3 kpc. This is likely an 
overestimate as the secondary should be H-poor, and the tem- 
perature of the viewed hemisphere is far from uniform, even 
at maximum. 

If we adopt the view that the X- and 7-ray emissions come 
from an energetic pulsar, we can make some additional esti- 
mates. Becker (2009) finds that pulsars have L x ~ 10~ 3 E, 
so from the observed X-ray flux we would infer a spin- 
down luminosity E pa (10 3 )47rd 2 2.4 x 10~ 13 ergcm _2 s _1 pa 
2.9 x 10 34 <i|_ c ergs _1 . Gamma-ray pulsars are also ob- 
served to follow a heuristic luminosity law L 7j fc e « ~ (E x 
lO^erg/s) 1 / 2 ( lAbdo et all 120101) . This is related to the 
observed flux as L 1 = Airfo J F 1 d 2 1 where the beaming- 
depen dent correction factor is /n ~ 0.7 — 1.3 ( Watte rTet al.l 
120091) . Thus the observed F 1 = 6.2 x ^"ergcm" 1 s _1 
gives a distance-dependent estimate for the spindown power 
of E « 5.5 x 10 34 /^c? 4 ergs _1 . This is consistent with the 
X-ray luminosity estimate when d pa 0.75/^ 1 kpc. However, 
if we match E to the ELC-estimated heating flux a distance 
of ~ 2 kpc is preferred. 

The nature of the intermittent 'flaring' emission is at present 
unclear. Since the bands were observed sequentially, we do 
not even know the characteristic color. All we can infer 
are large increases above the quiescent flux (up to ~ 5x 
in the optical; ~ 6x in the X-ray, reaching a peak flux 
~ 6 x 10 32 <i 2 pc ergs _1 ). These are concentrated toward the 
heated side of the companion, but visible at all phases. Shocks 
in the wind stripping the companion could be the origin. How- 
ever, since we do not know of such flaring activity on other 
well-studied black widow binaries, it is tempting to associate 
the activity wit h the unusually short Pb o f J 1311. As has 
been noted (eg. Richards & Albright 1993), low mass stars 
in close binaries often show enhanced coronospheric activity 
since tidal locking enforces rapid rotation, driving dynamo- 
generated magnetic fields. The other main ingredient for 
high magnetic activity is convection. This may suggest a low 
stellar-type mass for the companion. If sub-stellar, perhaps 
deep differential heating may drive the convection. A spec- 
troscopic search for coronal line emission, or spectra during 
an outburst, could resolve the origin. 

Additional data can greatly improve estimates of the sys- 
tem properties. Simultaneous multicolor photometry cover- 
ing many periods and flaring events can help us understand 
the temperatures of the various emitting regions. A complete 
spectroscopic study of the secondary can place important con- 
straints on the compact object mass and probe the state of the 
secondary photosphere. Spectra can also improve knowledge 
of the compact object's period, orbital epoch and semi-major 
axis, each of which are crucial for decreasing the parameter 



space to be searched for a 7-ray pulses. Of course, a 7-ray or 
radio pulsar ephemeris would produce a qualitative change in 
our knowledge of the system properties. We are pursuing all 
of these follow-up investigations. 
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Curvature Significance 

FIG. 4. — 2FGL bright sources in the spectral curvature - variability plane. 
Note the excellent pulsar(PSR)/blazar(AGN) separation. Binaries, a globular 
cluster (GC1), the Galactic Center (GC) and the unidentified sources belong 
to the pulsar group. 

We close with a comment on 7-ray pulsar population in the 
bright sample. Figure 4 shows these sources in the curvature- 
variability plane. The pulsar-blazar separation is excellent; 
the only pulsar in the blazar region is the Crab, whose PWN 
contributes the variable power-law component. The pulsar re- 
gion contains a few other source types: high mass binaries, 
a globular cluster, SNR, but each may harbor spin-powered 
pulsars. The unidentified sources, too, have significant cur- 
vature and low variability. Thus it is unsurprising that two 
of the six unidentified sources are good candidates for binary 
MSR although they lack radio detection. The other uniden- 
tified sources may be similar or new source types. An in- 
teresting question is whether J2339— 0533 and J 1 3 1 1 — 3429 
represent the first of a new source class - radio quiet MSR 
It is not yet clear that we need to appeal to such novelty: for 
both the optical gives evidence of a strong evaporative wind. 
It is entirely possible that radio emission is directed at Earth, 
but the pulse is scattered and dispersed to undetectability by 
this wind. Thus it remains possible that the MSP 7-ray beam 
lies within the radio emission zone - i.e. that all 7-ray MSP 
are radio detectable. 

Even if these two sources represent a harbinger of a new 
source class we emphasize that this source class must be 
small, at least for the bright po pulation, in contrast to some 
suggestions in the literature ([H arding, U sov & Mus limov 
12005b iTakata. Wang & Chengl 1201 ll) . The statistics in our 
sample clearly make this case - of the 41 young pulsars in 
this set only 17 (1/3) are radio selected. In contrast there are 
12 MSP in the bright sample, all radio pulsars. Thus even // 
all six unidentified sources prove to be radio-quiet MSP the 
radio loud subset will be at least 2/3 of the 7-ray MSP, twice 
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the fraction of the young pulsars. Thus, while J131 1—3429 
may have joined the 'black widow' club of pulsars evaporat- 
ing their companions, it is not yet clear if there is a 'radio quiet 
MSP' club for it to join. And if such a society does exist, it is 
very select, indeed. 
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